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ABSTRACT 

We have mapped the proto-binary source IRAS 16293-2422 in CO 2-1, ^^CQ 
2-1, and CO 3-2 with the Submillimeter Array (SMA). The maps with resolution 
of l'.'5-5" reveal a single small scale (~3000 AU) bipolar molecular outflow along 
the east-west direction. We found that the blueshifted emission of this small 
scale outflow mainly extends to the east and the redshifted emission to the west 
from the position of IRAS 16293A. A comparison with the morphology of the 
large scale outflows previously observed by single-dish telescopes at millimeter 
wavelengths suggests that the small scale outflow may be the inner part of the 
large scale (~15000 AU) E-W outflow. On the other hand, there is no clear 
counterpart of the large scale NE-SW outflow in our SMA maps. Comparing 
analytical models to the data suggests that the morphology and kinematics of 
the small scale outflow can be explained by a wide-angle wind with an inclination 
angle of ~30°-40° with respect to the plane of the sky. The high resolution CO 
maps show that there are two compact, bright spots in the blueshifted velocity 
range. An LVG analysis shows that the one located 1" to the east of source A 
is extremely dense, n(H2) ~10^ cm~^, and warm, Tkin > 55 K. The other one 
located 1" southeast of source B has a higher temperature of Tkm > 65 K but 
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slightly lower density of n(H2) ~10^ cm~^. It is likely that these bright spots are 
associated with the hot core-like emission observed toward IRAS 16293. Since 
both two bright spots are blueshifted from the systemic velocity and are offset 
from the protostellar positions, they are likely formed by shocks. 

Subject headings: stars: formation — ISM: individual (IRAS 16293-2422) — 
ISM: jets and outflows 



1. Introduction 



Outflow phenomena have been recognized as an important phase in the star forma- 
tion processes. It is generally accepted that when stars are formed by gravitational infall, 
outflo ws transfer excess angular momentum out of the system (e.g.. lShu. Adams, fc Lizano 
19871 ). Outflows are observed at various wavelengths from the radio to the ultraviolet, among 
which the millimeter /submillimeter bands probe the molecular outflows. Molecular outflows, 
recognized as high-velocity wings in CO and other molecular lines, are considered to be the 
ambient gas entrained or pushed by protostellar winds. By observing molecular outflows, we 
can determine their physical, kinematic and chemical properties, trace the mass-loss his tory 
of protostars, and probe the early phases of the star formation process (lArce et al.ll2007l . and 
references therein). 

IRAS 16293-2422 (hereafter referred to as 116293) is a well-studied proto-binary system 
located in the nearby molecular clou d complex L 1689. The distance to L 1689 is often 
assumed t o be 160 pc (IWhittetlll974l). howeyer som e recent studies suggested a distance of 



120 pc (Ide Geus et al 



1989 



Knude fc H0glll998l ). In this pap er we adopt 120 pc as the 
distance. The projected separation of the two protostars is ~5" (IMundy et al.lll992l ). which 



then corresponds to 600 AU. 116293 is classified as a Class young stellar object (YSO) since 
this system is undetectable at wavelengths short er than 10 Atm and shows a spectral energy 



distribution (SED) with very low temperature ( Andre &: Monterld Il9941). The bo lometric 
luminosity of 116293 is estimated to be ~36 Lp, (Correia. Griffin, fc Sara cenol 12004 ). 116293 



lummobiiy oi iiuzao ib esiimaieu lo ue ^ou uv^uiieia. vjiiimi. c<c oaiacenui^uu^JJi . i±uzyo 
is known to have a quadrupolar outflow ( Walker et al. 19881 : Mizuno et al. 199ol ): one bipolar 



pair extends along the east (blue) and west (red) (hereafter the El-W pair) directions, while 
the other one has an orientation from northeast ( red) to southwest ( blue) (h ereafter the 



NE-S W pair) with an inclination angle of 30°-45° ( iHirano et al.ll200ll ) or 65° ( IStark et al. 



20041 ). The NE-SW pair is more collimated, and its axis crosses the position close to the SE 
continuum source IRAS 16293A (hereafter source A). Therefore this pair has been tho ught 



to be powered by source A (e.g.. IWalker et al.lll988l : iMundy. Wootten. fc Wilkindll990l ). On 



the other hand, the E-W pair is not as well collimated and its axis, though not well defined. 
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seems to pass to the north of source A, where the continuum source IRAS 16293B (hereafter 
source B) is the only possible dri ving; YSO. Hence source B has often been assumed as the 
driving source of the E-W pair. IStark et al.l (120041 ) showed that source B has very narrow 
line widths, a low luminosity, and is not evidently associated with high-velocity gas. They 
interpreted source B as a T-Tauri star that drove the E-W pair and that the E-W pair is 
now a fossil flow. However, these interpretations are based on observations at low angular 
resolution (i.e., >5"), that are insufficient to resolve the two protostars and reveal the relation 
between the quadrupolar outflow and the binary. 



In fact, using higher angular resolution observations IChandler et al.l (120051 ) proposed an 
alternative interpretation. They resolved source A into four components, two at 1 mm and 
another two at 7 mm, and suggested that two o f them are responsible for driving the NE-SW 
and the E-W pair, respectively. In addition, IChandler et al. (12005^ proposed that source 
B is actually much younger than previously interpreted and may have not yet begun the 
phase of mass loss. Their proposition was, however, made based on (1) the dust continuum 
morphology and proper motions of source A, and (2) the velocity structures of emission of 
H2CO, H2S, and SO, rather than CO, which is a much better tracer of the overall outflow 
structure free from chemical peculiarities and shocks. 

In order to study the molecular outflows in the vicinity of 116293 in more detail, we 
have carried out CO 2-1, ^'^CO 2-1, and CO 3-2 observations at higher angular resolutions 
using the Submilhmeter Arrajfl (SMA). Our SMA observations provide sufficient angular 
and spectral resolution in the CO lines to resolve the kinematics and the structure of the 
outflow in the central region of 116293. In this paper, we summarize the observational details 
in §2; present our SMA CO 2-1, ^^CO 2-1, and CO 3-2 results in §3; discuss the outflow 
driving mechanism, physical parameters of the outflow, and the kinematics in the vicinity of 
116293 in U. 



2. Observations 



The Submilhmeter Array (SMA) (IHo. Moran. fc Ldl2004l ) observations of the CO J=2-l 
(230.538 GHz), CO 3-2 (345.796 GHz), and ^^CO 2-1 (220.399 GHz) transitions were carried 
out during February 2003 to February 2005. The phase tracking center was q;(J2000) = 16''32™22.91' 
(5(J2000) = -24°28'35.52". The bandwidth of the spectral correlator was 656 MHz in 2003 



^The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and 
the Academia Sinica. 
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and 2 GHz in 2005. The frequency resolution used was 812.5 kHz, corresponding to the 
velocity resolution of 1.06 km s~^ at 230 GHz and 0.70 km s~^ at 345 GHz. The field of view 
of the array is ~54" at 230 GHz and ~36" at 345 GHz. The CO 2-1 and ^^CO 2-1 lines 
were observed simultaneously on 18 February 2005 with 6 antennas, providing 15 baselines 
ranging from 10 to 70 m. The CO 3-2 data were observed on February 26 and April 23 in 
2003 with 4 and 5 antennas, respectively, providing 10 independent baselines ranging from 
12 to 120 m. We note that the CO 3-2 observations provide fewer short baselines than the 
CO 2-1 observations, therefore the CO 3-2 data are less sensitive to extended structures. 

The visibility data were calibrated using the MIR softw are packa ge originally developed 



for Owens Valley Radio Observatory (IScoville et al.l Il993l : lOi 



20051 ) and modified for the 



SMA. The absolute flux calibration and bandpass calibration were done by observing Callisto 
(on 2003 February 26), GanjTiiede (on 2003 April 23), and Uranus (on 2005 February 18). 
The uncertainty in the absolute flux scale is estimated to be ~10%. The visibility phase 
and amplitude calibration was achieved by observing the quasars NRAO 530 (J1733-130) 
in 2003 and J1743-038 in 2005. The continuum level was measured from line-free channels 
in the calibrated visibility data. The MIRIAD and the NRAO AIPS packages were used to 
make the spectral images and for further analysis. The maps were restored using uniform 
weighting, resulting in synthesized beam sizes of 4'.'7x2'.'5 with a position angle of 0.0° for 
the CO 2-1, 4'.'9x2'.'6 with a position angle of 1.5° for the ^^CO 2-1, and l'.'6xl'.'4 with a 
position angle of 26.2° for the CO 3-2. The rms noise levels of the channel maps with ~1 
km s~^ width are ~0.15 Jy beam^^ for the CO 2-1, ~0.13 Jy beam~^ for the ^^CO 2-1, and 
~0.3 Jy beam-i for the CO 3-2. 



3. Results 
3.1. CO 2-1 

Fig. [1] shows the velocity channel maps of the CO 2-1 emission uncorrected for the 
primary beam attenuation. The CO 2-1 emission is shown in the channels from Vlsr = 
—7.2 km s~^ to 31.9 km s~^. In most of the channels, the spatial distribution of the CO 2-1 
emission is elongated along the east-west direction. It is likely that this elongat ed CO emis- 



sion originates from the molecular outflow previously imaged on larger scales (jStark et al. 



20041 ). The CO emission at Vlsr = 4.4 km s ^ is scattered over the field of vi ew. Since this 



veloc ity channel is close to the cloud systemic velocity of Vlsr = 4.0 km s ^ (IMizuno et al. 



19901 ). most of the emission has been resolved out due to the lack of short baselines. The 
CO distribution at Vlsr = 3.3 km s~^ and Vlsr = 5.5 km s~^ shows the east-west elon- 
gation, suggesting that the CO emission at these channels comes from the outflowing gas 
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even though the velocity offset from the systemic velocity is smalL Therefore, we infer that 
the CO 2-1 emission seen in the channels from Vlsr = —7.2 km s~^ to 3.3 km s~^ is the 
blueshifted outflow component and that in the channels from Vlsr = 5.5 km s^^ to 31.9 km 
s~^ is the redshifted outflow component. 

The blueshifted emission is extended to the east of source A and the west of source B as 
well. The eastern component is seen in most of the blueshifted channels from Vlsr = —7.2 
km s~^ to 3.3 km s~^, while the western component is in the velocity range from Vlsr = 
—0.9 km s^^ to 3.3 km s~^, and is closer to the systemic velocity compared to the eastern 
component. As shown in the velocity channels from Vlsr = —0.9 km s~^ to 1.2 km s~^, the 
eastern blueshifted component has a fan shape pointing to source A with a position angle of 
-75°. 

The redshifted emission extends to the east and the west of source A. The eastern 
redshifted component mostly overlaps with the eastern blue emission. The western redshifted 
component has a fan shape structure with its apex at the position of source A at P. A. ~90°. 
This structure appears from Vlsr = 8-6 km s~^, and the emission moves away from the 
source as the velocity increases. 

The total integrated intensity (zeroth moment) maps of the blueshifted and redshifted 
CO 2-1 line emission are presented in Fig. [2^ (uncorrected for primary beam attenuation). 
In the blueshifted component, there are two prominent peaks 1" to the east of source A 
(labeled bl) and 1" southeast of source B (labeled b2). The peak brightness temperature of 
CO 2-1 emission at bl and b2 are 57 K and 63 K, respectively. These two peaks are seen 
in all the blueshifted channels (Fig. [1]) starting from Vlsr = —6.2 km s~^. We note that 
both bl and b2 have not been resolved in the previous CO outflow maps observed with the 



single-dish telescopes ( IStark et al.ll2004j ). 



In Fig. [3^ we compare the CO 2-1 spectrum observed by the SMA with that observed 
by the JCMT. We have applied the primary beam correction, smoothed the SMA map to 
have a 21" Gaussian beam, which is same as the JCMT beam at 230 GHz, and obtained 
the spectrum at the phase tracking center. The JCMT spectrum was resampled to have the 
same velocity resolution as the SMA data. To convert the antenna temperature Ta* into 
brightness temperature, we adopted a main beam efficiency, ?7mb, of 0.69 for the JCMT at 
230 GHz. It is shown that the SMA recovered most of the flux in the velocity range lower 
than Vlsr ~2 km s~^, indicating that most of the blueshifted emission in this velocity range 
comes from the compact regions labeled bl and b2. The SMA spectrum is slightly brighter 
than the JCMT spectrum in several channels, probably due to the calibration uncertainties 
in both the SMA and JCMT data. In the redshifted velocity range, the SMA recovered 
~70% of the flux at Vlsr > 10 km s~^. The flux in the velocity range close to the systemic 
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velocity is recovered ~60% except the channel of the deep absorption dip at Vlsr = 4.4 km 
s^^ At this velocity, absorption caused by the foreground cloud also contributes to the total 
absence of emission in the interferometric map because the JCMT spectrum also shows a 
deep absorption dip at the same velocity. 

3.2. CO 3-2 

The CO 3-2 total integrated intensity (zeroth moment) map is shown in Fig. [2]d (uncor- 
rected for primary beam attenuation). The emission is integrated over the velocity channels 
from Vlsr = —7.5 km s~^ to 3.5 km s~^ on the blueshifted side, and from Vlsr = 5.5 km 
s~^ to 31.5 km s~^ on the redshifted side. There is no significant emission in the velocity 
channel centered at Vlsr = 4.5 km s~^. The CO 3-2 emission in this velocity range arises 
from the extended ambient cloud and is poorly sampled by the uv coverage of the CO 3-2 
data. There are two bright blueshifted components, one located 2" southeast of source A 
and the other located 1" southeast of source B, that are likely to be the counterparts of bl 
and b2 identified in the CO 2-1 map (see Fig. [2^). The peak brightness temperature of 
CO 3-2 emission is 39 K at bl and 55 K at b2. These two components are connected by 
faint emission. The blueshifted emission extending toward the northeast with a P.A. of ~70° 
corresponds to the eastern blueshifted emission in the CO 2-1 map. As in the case of CO 
2-1, the redshifted emission extends to the east and west of source A. 

In order to estimate the missing flux, we have smoothed the SMA CO 3-2 map to 
14" resolution, which is the HPBW of the JCMT at 345 GHz, and compared the spectrum 
at the phase tracking center with that obtained by the JCMT (Fig. ^jp). The primary 
beam correction was applied to the SMA CO 3-2 map before it was smoothed. The JCMT 
spectrum was converted into the brightness temperature unit by assuming t^mb = 0.63 at 
345 GHz. Fig. [3]d shows that the missing flux of the SMA CO 3-2 data is larger than that 
of the CO 2-1 except for the high- velocity part of the redshifted wing at Vlsr ^ 14 km s~^, 
where most of the flux is recovered by the SMA. The SMA recovered ~50% of the CO 3-2 
flux in the velocity ranges —9.0 < Vlsr < 2.0 km s"^ and 6.0 < Vlsr < 13.0 km s~^. That 
more missing flux is evident in the CO 3-2 data than the CO 2-1 data is likely due to it 
containing fewer short baselines, as mentioned in §2. 
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3.3. ^^CO 2-1 

Fig. [3t compares the ^'^CO spectrum observed with the SMA and that observed with 
the JCMT. The SMA spectrum was obtained after we corrected for the primary beam at- 
tenuation and smoothed to 21" resolution. Due to a software problem, the velocity values 
of the ^^CO 2-1 data were not recorded correctly. We found that the ^^CO 2-1 spectrum 
observed with the SMA was blueshifted by ~1 km s~^ with respect to that observed with 
the JCMT. Therefore, we calibrated the velocity of each channel in the SMA data using 
the ^^CO 2-1 spectrum observed with the JCMT. Most of the ^^CO 2-1 flux is recovered 
by the SMA at Vlsr < 2.0 km s~^. The velocity channel maps of the ^^CO 2-1 emission 
presented in Fig. H] (uncorrected for primary beam attenuation) show that emission in this 
velocity range mainly comes from the compact source close to source A. It is likely that 
this component corresponds to bl in the CO 2-1 map. On the other hand, the brightest 
component b2 in the CO 2-1 map is less significant in the ^^CO 2-1. It is identified as a 
northwestern extension of the emission in the velocity channels from —0.8 to 1.5 km s~^. The 
^^CO emission extends to the east of source A in the redshifted velocity channels at Vi^sr 
= 5.9 km s^^ and 7.0 km s~^, and extends to the east and west of source A at Vlsr = 8.1 
km s~^. It should be noted that a similar structure is seen in the velocity channel maps of 
the CO 2-1 (see Fig. [1]), suggesting that the ^^CO 2-1 emission in this velocity range arises 
from the dense part of the outflow. The ^^CO 2-1 flux recovered by the SMA is 60-70% 
in the velocity channels from Vlsr = 5.9 km s^^ to 9.2 km s~^. In the velocity channels at 
Vlsr = 2.6 km s~^ and 3.7 km s~\ the ^^CO 2-1 map shows a V-shaped structure. Since 
the southwestern and southeastern edges of the ^^CO 2-1 correspond to the northern edges 
of the western and eastern lobes, respectively, it is likely that the ^^CO 2-1 emission in this 
velocity range comes from the envelope excavated by the outflow. Note that the primary 
beam correction is not applied to the images shown in Fig. HI 



4. Discussion 
4.1. Outflows in the vicinity of IRAS 16293-2422 

4.1.1. Driving source 

As was shown in §3, the CO 2-1 and 3-2 maps observed with the SMA clearly show 
small scale outflow emission in the vicinity of 116293. It is quite natural to assume that the 
CO 2-1 and 3-2 outflow emission has the same origin. Each of the CO 2-1 and 3-2 outflows 
has a blueshifted component extending to the east and a redshifted component extending to 
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the west of source A (see Fig. [2^ and[2)D). In the CO 2-1 channel maps, both the eastern 
blueshifted component and the western redshifted component show fan-shaped patterns with 
their apexes pointing toward source A. Ahhough the eastern redshifted component is less 
clearly extending from the position of source A, its coincidence with the eastern blueshifted 
emission suggests that the eastern blue- and redshifted components represent the front side 
and the back side of a single outflow lobe, respectively. Such an outflow lobe showing both 
blue- and redshifted emissions together is li kely to be an outflow wi th a wide opening angle 



and an axis close to the plane of the sky (jCabrit fc BertoutI Il986l ) . These morphological 



characteristics of the small scale outflow suggest that it is most likely driven by source A. 

The direction of the outflow axis is rather difficult to define clearly. The western red- 
shifted lobe in CO 2-1 shows a position angle of almost —90°, while the blue- and redshifted 
eastern lobes show position angles of ~75°. The axis of the CO 3-2 outfiow is approximately 
along a position angle of 80°. These results suggest that the axis of the small scale outfiow 
is approximately the east-west direction. The difference in direction between the blue- and 
redshifted components of the CO 2-1 emission, and the CO 2-1 and 3-2 emission, could also 
indicate that the high velocity gas comprises emission from two or more outfiows. There 
may be an interaction between the wind or jet driving the large scale NE-SW outfiow and 
the molecular gas, as indicated by the offset of shock-tracing molecule line emission peaks 
to the NE of source A (Chandler et al. 2005). However, the bulk of the CO emission we 
observe is consistent with a predominantly E-W outflow. 



4-. 1.2. Comparisons with the large scale outflows 



It is important to know the relationship between the small scale (~3000 AU) outflow 
observed with the SMA and the large scale (~15000 AU) outflows previously observed with 
single-dish telescopes even though the size scales are quite different between them. We should 
note that the missing flux and the single fleld of view in the SMA observations prevent us 
from directly relating the two types of outflows on these very different scales. Nevertheless, 
morphological comparisons of the small and large scale outflows could allow us to infer their 
relationship. 

Single-dish observati ons showed two large scale outflows; one is the E-W pair and the 
other is the NE-SW pair JWalker et al.lll988l : ]Mizuno et al.lllQQol : Ihis et al.l[2002[lstark et al 
20041 ). The blue- and redshifted lobes of the E-W pair are located to the east and the west, 
respectively, showing the same velocity trend as the small scale outflow. The outflow axes 
of these two outflows are also similar. These facts naturally suggest that the small scale 
outflow is the inner part of the large scale E-W pair, and that the bulk of the CO emission 
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observed by the SMA is due to the E-W outflow. If this is the case, then the E-W pair is 
not a fossil flow but a currently active one, and the driving source of the E-W pair is most 
likely to be source A. 



Do the SMA maps show any counterparts of the NE-SW pair? IChandler et al.l (120051 ) 
proposed that the driving source of the NE-SW pair is one of the radio sources, A2, based 
on its bipolar morphology along the NE-SW direction in their 43 GHz map. However, our 
SMA maps show no counterpart of the NE-SW pair, if we consider the small scale outflow 
mapped with the SMA as the inner part of the E-W pair. One possible reason for why there 
is no counterpart of the NE-SW pair in the SMA maps is that the NE-SW outflow lobes are 
l ocated outside the SMA field of vi e w at both 230 an d 345 GHz. In fact, single-dish maps 
dMizuno et aDll990l : Ihis et al.l[2002l : Istark et al.1l2004l ) show that both blue- and redshifted 
lobes are detached by ~1' from the centroid of 116293. This implies that the NE-SW pair is 
older than the E-W pair, so the NE-SW pair has already cleaned out the gas along its travel 
direction. However this interpretation is totally different from those of previous studies, 
where the NE-SW pair is considered younger than the E-W pair. On the other hand, the 
NE-SW pair is better collimated than the E-W pair. Observations of CO outflows of YSOs 
at different evolutionary stages suggest that the opening angl e of the outflow lobe becomes 



larger as the central source evolves (e.g.. lArce fc Sargentll2006l ). If we take this into account 



it is natural to consider that the well-coUimated NE-SW pair is younger than the E-W pair. 
It is possible that the small scale counterpart of the NE-SW pair is absent because this 
outflow is currently inactive due to its episodic nature. E posodic outflows have been found 



in several YSOs, such as L1157 (e.g., iGueth et al.lll996l ). Since there is no counterpart of 



the NE-SW pair in our SMA maps, we are unable to determine if the NE-SW outflow is 
driven by source A of B, and to support or refute the diff erent claims by previous studies 
regarding this flow (jStark et al.ll2004l : IChandler et al.ll2005l ). 



4.2. Kinematics of the small scale outflow 



Many models have been proposed to explain the relation between molecular outflows 
and protostellar jets. Of these, two scenarios are considered plausible: (1) the jet-driven 
model, w here molecular outflows consist of ambierit gas swept up by the bow shock of the 
jet head jRaea fc Cabritlll993l : Masson fc CherninI Il992l : Ichernin et al1ll994h . and (2) the 
wind-driven mode L where molecular outflows are the swe pt-up gas entrained by wide-angle 
magnetized wind (IShu et al.l Il99ll . l2000l : iLi fc Shul Il996l ). These scenarios have different 
kinematical features in a Position-Velocity (P-V) diagram obtained from a cut along the 
outflow axis. The jet-driven model will produce a convex structure in such a P-V diagram 
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with high velocity components at the jet head; while the wind- driven model w ill introduce 
a parabolic structure that originates from the central star (e.g. 



LeeetaLl 120001). 



In a CO 2-1 P-V diagram, extracted along the outflow axis at P. A. = 90° (Fig. 
we found that the western redshifted lobe shows a parabolic shape, as is expected from the 
wind-driven mode l. We therefore a dopted the simplified analytical model of a wind-driven 
shell proposed by iLee et al.l (l2000l ) to study the kinematics of the small scale outflow. In 
the cylindrical coordinate system, the structure and velocity of the shell can be written as 
follows: 



Z = CR , Vr = VqR, Vz = VqZ, (1) 

where z is the distance along the outflow axis; R is the radial size of the outflow perpendicular 
to z ; C and Vq are free parameters that describe the spatial and velocity distributions of 
the outflow shell, respectively. 



First, we generated the model curve that dehneates the outflow shell featurcl (see Fig. [6]), 
to determine the free parameter C. In Fig. [H] the green line and the black line delineate the 
curves produced by C = 0.20 arcsec"^ and C = 0.10 arcsec"^, respectively. Both curves 
assume the inclination angle of the axis from the plane of the sky, i to be 30°. However, 
the curve is insensitive to the inclination angle if i is smaller than 40°. Since both eastern 
and western lobes have opening angles of 60°-80°, and are accompanied by the CO emission 
component with opposite velocity shifts (i.e. the eastern blue lobe is accompanied by the 
redshifted component and the western red lobe by the blueshifted component), it is unlikely 
that the inclination angle is larger than 40°. Next, we determined vq and i by comparing 
the velocity pattern of the model curve with the observed P-V diagram. The green curve in 
Fig. El was produced hj i = 30°, C = 0.20 arcsec"^, and vq = 1.3 km s"^ arcsec"^; while the 
dashed curve was produced by i = 40°, C = 0.10 arcsec"^, and Vq = 0.8 km s~^ arcsec"^. A 
wind-driven outflow model predicts the "Hubble law" velocity str ucture, that is the velocity 



incre ases as the emission moves away from the central source (e.g.. lShu et al.lll99ll : iLee et al. 



2000l ). As we have mentioned in §3.1, the redshifted lobe of the small scale outflow has this 
signature. Hence, our simple modeling qualitatively indicates that the driving mechanism 
for the small scale outflow could be a wide-angle wind inclined to the plane of the sky by 
~30° to 40°. The dynamical age of the small scale outflow is given by 1/vo and is 450-700 
years. 

It should be noted that the wind-driven model does not contradict to the evidence 



^We only model the red lobe because it shows a clear fan-shape morphology. 
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of a precessing jet presented by I Chandler et al.l (120051 ) . Recent numerical simulations of 
magnetocentrifugal winds pr edict an on-axis density enhanceraent wi thin the X-wind type 
of wide opening angle wind (ILi &: Shul Il996l : IShang et al.l |2002| . |2006| ) , and so it is possible 
to have a jet-like feature in such a model that does not show up in CO observations. 



4.3. Physical parameters of the small scale outflow 

We estimated the outflow mass using the CO 2-1 and ^'^CO 2-1 data cubes corrected 
for the primary beam attenuation and assuming LTE conditions. Since the peak brightness 
temperature of the CO emission reaches ~50 K, we assumed 50 K as an excitation tempera- 
ture. We adopted an H2 to CO abundance ratio of 10^ and a mean atomic weight of the gas 
of 1.36. A comparison between the channel maps of CO 2-1 (Fig. [1]) and those of ^'^CO 2-1 
(Fig. H]) suggests that the CO 2-1 emission is significantly optically thick in the low- velocity 
redshifted component at Vlsr ~5-9 km s~^ and the low-velocity blueshifted component at 
Vlsr ~0-2 km s^^. The ^^CO 2-1 emission is also detected in the positions of the bright 
spots bl and b2; however, we excluded these components from our calculation. The physical 
conditions of these bright spots are discussed in the next subsection. At velocities where the 
^^CO emission is detectable, we estimated the optical depth of the CO 2-1 line using the 
following equation: 

^^(CO) [CO] /1-e 



whe re [CO]/[^'^CO] is the abundance ratio between CO and ^^CO that is assumed to be 



77 (IWilson &: Roodlll994j ). The optical depth of the CO 2-1 emission in the low- velocity 
redshifted component is estimated to be ~3-36 in the western lobe and ~6-24 in the east- 
ern lobe, and that of the low-velocity blueshifted component is ~2-17 in the eastern lobe. 
Outflow masses in the optically thick velocity channels were corrected using: 

Mthick = ^^mivi) = m'{vi) (^y^^T^i) ' 

where m'{vi) is the mass at velocity Vi without optical depth correction; m(i;j) is the mass 
after optical depth correction; r, represents the optical depth of CO 2-1 at velocity Vi. 
Outside the velocity ranges where the ^^CO emission was observed, we assumed that the 
CO emission is optically thin, and obtained the mass by Mthm = X] ^'{'^i)- The total mass 
of the outflow was obtained by Mtotai = ^thin + Mthick- Table [T] summarizes the masses of 
the blueshifted and redshifted components in the eastern and western lobes. The optically 
thick part of the spectrum dominates the outflow mass calculations. The outflow masses are 
underestimated by an order of magnitude without the correction for optical depth. Since 
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significant flux is missed in the low velocity range due to the lack of short baseline data 
and absorption by foreground clouds, the masses derived here are lower limits. In both 
eastern and western lobes, the redshifted component is signiflcantly more massive than the 
blueshifted component. This suggests that dense ambient gas exists behind the outflow 
lobes, so that the blueshifted component in the front can expand rather freely, while the 
redshifted component in the back needs to push more material. 

The dynamical parameters of the flow are summarized in Table 2. Table 2 gives two 
kinds of values; one is corrected for the inclination effect {i = 35°) and the other is without 
this correction (uncorrected). The inclination angle of 35° is chosen based on the kine- 
matical modeling described in the previous section. In order to compare the energies of 
the small scale outflow with those of the large scale E-W outflow observed with single-dish 
telescopes, we estimated the momentum supply rate and the stellar mass loss rate. These 
dynamical quantities can be derived reasonably well even though the map is incomplete (e.g. 



Bontemps et al.lll996l ). The momentum supply rate and the mass loss rate of the small scale 



outflow lobes are estimated to be ~2xl0 Mq km s ^ yr ^ a nd 9.0x10 M^t^ yr ^, respec- 



tively. The CO 1-0 single-dish observational result of 116293 (iMizuno et al.lll990l ) scaled to 
a distance of 120 pc suggests that the momentum supply rate and the mass loss rate of the 
large scale E-W outflow are ~1.6xl0~^ Mq km s~^ yr~^ and ~8.0x 10"''' Mq yr~^ in the 
eastern lobe, while those of the western lobe are ~l.lxlO~^ Mq km s^^ yr^^ and ~5.6xl0~'' 
Mq yr~^, respectively. Although the parameters we derived from the SMA data are lower 
limits, it is likely that the small scale outflow has comparable momentum and energy with 
the large scale E-W outflow. This may imply that the E-W outflow has been fairly constant 
in time, if the compact outflow is the inner part of the large scale E-W outflow. 



4.4. Nature of the bright, compact components 



As shown in §3, the high resolution CO maps reveal two bright and compact com- 
ponents, labeled bl and b2, at blueshifted velocities (see Fig. [T] and F ig. [2]). In order to 
inves t igate the nature of bl and b2, we carried out an LVG c alculation (IGoldreich fc Kwan 
1974 : IScoville fc Solomonlll974l ) using the code developed by lOka et al.l (119981 ). The veloc- 



ity gradient dv/dr is assumed to be 1000 km s~^ pc~^ as estimated from the CO 2-1 line 
width of ~ 7 km s~^ and the size of ~0.005 pc of bl and b2. The LVG calculation result is 
shown in Fig. [71 Fig. [7^ shows the iso-intensity contours of the CO 2-1 line (solid contours) 
and iso-line ratio contours of CO 2-1/^^CO 2-1 (dashed contours) as a function of kinetic 
temperature and H2 number density. Fig. [Tb shows the CO 2-1 iso-intensity contours (solid 
contours) and iso-line ratio contours of CO 2-1/^'^CO 2-1 (dashed contours). 
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In order to compare the CO 3-2, CO 2-1, and ^^CO 2-1 lines at the positions bl and 
b2, we made the maps of these three hnes using the visibihty data in the range from 10 kA 
to 52 kA, corresponding to the uv range common for the three data sets. The beam sizes in 
uniform weighting of the CO 3-2, CO 2-1, and ^^CO 2-1 maps were 3'/lx2'.'6 (P.A. 5.2°), 
4'.'7x2'.'5 (P.A. -3.0°), and 4'.'8x2'.'6 (P.A. 1.5°), respectively. In spite of the same uv range 
limit, the beam size of the CO 3-2 map is different from those of the CO 2-1 and ^^CO 
2-1 due to differet uv coverages. Therefore, we convolved the three maps to the same beam 
size, i.e. 4'.'9x2'.'7 with a P.A. of 0.0° after the correction of the primary beam attenuation. 
The CO 3-2, 2-1 and ^^CO 2-1 spectra obtained at the positions of bl and b2 are shown 
in Fig. [HI The brightness temperature of CO 2-1 at the position of bl is Tb ~50 K, and 
that at b2 is Tb ~60 K. The CO 2-1 flux integrated over the velocity range from Vlsr = 
-8 km s-^ to 2 km s^^ is 102.2 K km s'^ at bl and 207.5 K km s'^ at b2. The ^^CO 2-1 
flux integrated over the same velocity range is 55.9 K km s~^ at bl and 18.7 K km s~^ at 
b2. Since most of the CO and ^^CO 2-1 fluxes in this velocity range were recovered by the 
SMA, the effect of missing flux should not be signiflcant. The CO 2-1/^^CO 2-1 intensity 
ratio of 1.8 at bl suggests that this component is warm (T^in ~55 K) and extremely dense 
(n(H2) ~10'^ cm-3). On the other hand, the CO 2-1/^^00 2-1 ratio at the position of b2 is 
~11, suggesting that b2 is slightly warmer (Tkm ~65 K) but less denser (?t-(H2) ~10^ cm~^) 
than bl. 

The CO 3-2 flux integrated over the same velocity range as CO 2-1 is 58.4 K km s~^ 
at bl and 105.6 K km s~^ at b2. The CO 3-2/2-1 ratio at the positions of bl and b2 
calculated from the integrated intensity values are 0.57 and 0.51, respectively. On the other 
hand. Fig. [7] shows that the line ratio should be >0.9 if the CO 2-1 intensity is brighter 
than 20 K. However, in contrast to CO 2-1, only ~50% of the CO 3-2 flux in this velocity 
range was recovered by the SMA. If we assume that most of the CO 3-2 flux within the 
central 14" region arises from the compact sources bl and b2, their actual flux values might 
be twice the calculated values. If we take this effect into account, the CO 3-2/2-1 ratio 
might be close to 1 at both bl and b2. Therefore, the CO 3-2 results do not contradict 
the physical conditions inferred for bl and b2 from the CO and ^'^CO 2-1 data. We should 
note that there are two effects that could contribute to a low CO 3-2/2-1 ratio. If the CO 
emission is assumed to be in LTE and optically thick, the CO 3-2 emission has a higher 
optical depth than CO 2-1. Therefore if there are temperature gradients due to internal 
heating sources, the CO 3-2 emission will tend to trace outer and cooler gas. Even for a 
uniform temperature, the r = 1 surface of the CO 3-2 emission with higher optical depth 
will be more extended and more flux will be resolved out by the interferometer, independent 
of the other uv coverage issues. 

It should also be noted that the kinetic temperature and density derived here are the 
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averaged values over the beam area of 4'/9x2'.'7. Since the bright spots bl and b2 are compact, 
the CO 2-1 brightness temperature at bl and b2 might be higher if we observe with higher 
angular resolution. Indeed, the brightness temperatures of CO 2-1 and CO 3-2 at bl and 
b2 measured in the higher resolution maps are higher than those of the spectra shown in 
Fig. [HI Therefore it is possible that the kinetic temperatures at bl and b2 are higher than 
those estimated here. 

Complex molecules, such as HCOOCH .s and C^H.^CN, have b een detected toward both 
source A and source B (IKuan et al.l |2004J : IChandler et al.l l2005l ). The estimated kinetic 
temperatures (>55 K at bl and >65 K at b2) are consistent with t he rotational temp erature 
of ~100 K derived from the analysis of the SMA HCOOCH3 data JKuan et al.ll2004h . Since 
the bright components bl and b2 are close to source A and B, respectively, it is likely that 
these warm and dense components are associated with hot core-like emission. Our CO 2-1 
and 3-2 maps show that both bl and b2 are blueshifted from the systemic velocity and are 
offset from the positions of continuum peaks. These results support the i dea that shocks play 



an important role in formin g hot cores around low mass protos t ars (e .g. Ivan Dishoeck et al. 



I995I : IChandler et al.l 120051 ) . As discussed in IChandler et al.l (l2005l ) , it is likely that an 
interaction between the jet from source A and extremely dense (~10^ cm~^) gas at bl 
contributes to the observed rich chemistry in this region. On the other hand, it is unclear 
how the b2 component is formed. This compone nt is also seen in HC N 4-3 emission at 
blueshifted velocities from —2.62 to +0.65 km s~^ fiTakakuwa et al.ll2007l). I n addition, high 



Chandler et al 



((20051) show the local 



resolution H2CO, SO, and H2S maps presented by 
emission peaks close to b2 rather than the continuum peak of source B. Since the sulphur - 
bearing molecules and H2CO are known to be enhanced in shocks (e.g. iBachiller et al.ll200ll ). 
the component b2 is also likely to be a s hocked re g ion. If sou rce B is a proto s tar in the 
earlier stage than source A as argued by IWoottenI (119891 ) and IChandler et al.l (120051 ) . b2 
may be a spatially compact outflow from source B that can interact with dense envelope 
gas and produce shocks. Another possibility is that b2 indicates the shock produced by the 
interaction between the outflow from source A and the envelope surrounding source B. 



5. Conclusions 

We have carried out CO 2-1 ^^CO 2-1, and CO 3-2 observations of the proto-binary 
system IRAS 16293-2422 with a resolution of l'.'5-5" using the Submillimeter Array (SMA), 
which are about a factor of 5 better than previous observations in these lines. Our obser- 
vations reveal the detailed structures of molecular outflows close to the binary system. Our 
main results are summarized as follows: 
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1. The high resolution images of CO 2-1 and CO 3-2 show a compact bipolar outflow- 
along the east-west direction in the vicinity of the binary. The center of this small scale 
outflow is close to source A, suggesting that source A is the driving source. If we interpret 
this small scale outflow as the inner part of the large scale E-W pair of the quadrupolar 
outflow, the E-W pair is a currently active outflow rather than a fossil flow. On the other 
hand, there is no clear counterpart of the large scale NE-SW pair in our interferometric 
maps. 

2. The shape and velocity structure of the western redshifted lobe arc well described by 
an analytical model suggesting that the small scale E-W pair could be driven by a wide-angle 
wind with an inclination angle of ~30°-40°. 

3. There are two compact and bright spots at blueshifted velocities. One is located 1" 
east of source A (labeled bl) and the other is 1" southeast of source B (labeled b2). An 
LVG analysis shows that bl is extremely dense, n(H2) ~10'' cm~^ and warm Tkin > 55 K, 
while component b2 has a higher temperature of Tkin > 65 K but slightly lower density of 
n(H2) ~10^ cm"^. It is hkely that these bright spots are produced by means of shocks and 
are associated with hot core-like activity observed toward both source A and source B. 

We would like to thank all the SMA staff who have helped us to carry out the observa- 
tions. We also wish to thank the anonymous referee who gave the insightful comments. N. 
H. is supported in part by NSC grant 96-21 12-M-001-023. Dr. Y. -N. Su is acknowledged 
for fruitful discussions on this research. 
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Table 1. Masses of the small-scale outflow 





Eastern lobe 


Western lobe 


Component 


(10-=^ Me) 


(10-=^ Mo) 


Blue'" 


5.4 


1.1 


Red"^ 


19.1 


12.5 


Total 


24.5 


13.6 



'"Velocity range: Vlsr= —7.2-3.3 km s ^ 
'^Velocity range: Vlsr= 5.5-31.9 km s~^ 



Table 2. Dynamical parameters of the small-scale outflow 



Eastern lobe Western lobe 



Parameters 


i=35° 


uncorrected 


z=35° 


uncorrected 


Momentum (M© km s""^) 


10.8x10-2 


6.2x10-2 


9.1 xlO-2 


5.2x10-2 


Kinetic energy (erg) 


12.1x10^2 


4.0x10^2 


21.9x10^2 


7.0x10^2 


Dynamical timescale (yr) 


600 


850 


500 


700 


Mass loss rate*^ (M0 yr""*^) 


9.0x10-^ 


3.6x10-^ 


9.0x10-^ 


3.6x10-^ 


Momentum supply rate (Mq km s~^ yr~^) 


1.8x10-^ 


7.2x10-^ 


1.8x10-^ 


7.3x10-5 


Mechanical luminosity (L©) 


0.31 


7.2x10-2 


0.42 


9.6x10-2 



^The central wind speed is assumed to be 200 km s ^. Then we derived the mass loss rate by 
assuming momentum conservation of the entrained CO gas and the central wind. 
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R. A. offset (arcsec) 

Fig. 1. — Velocity channel maps of the CO 2-1 emission. Offsets (in arcseconds) are with 
respect to the phase tracking center, a = 16''32™22.91^ and 6 = -24°28'35.52"(J2000). The 
contour levels correspond to —5a, 5a, 10a, 20a, then increase in every 20a step. The la 
level is 0.15 Jy beam^^. The crosses indicate the positions of the continuum sources A and 
B; the open circle is the SMA field of view at 230 GHz, and the filled ellipse indicates the 
synthesized beam. Note that the primary beam correction is not applied to the data. 
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Fig. 2. — a: Integrated blueshifted (blue contours) and redshifted (red contours) CO 2-1 
maps. The emission is integrated over the velocity channels from Vlsr = —7.2 km s^^ to 3.3 
km s~^ on the blueshifted side and from = 5.5 km s^^ to 31.9 km s~^ on the redshifted 
side. Contours are drawn at —5a, 5a, 10a, 20a, then every 20a step. The la levels are 0.6 Jy 
beam~^ km s'^ for the blueshifted component and 0.9 Jy beam"^ km s~^ for the redshifted 
component. The two crosses denote the positions of the two continuum sources. Positions of 
the two prominent compact components are labeled (bl and b2). The open circle indicates 
the field of view at 230 GHz, and the filled ellipse indicates the synthesized beam size, b: 
Same as Fig. 2a for the CO 3-2. The emission is integrated over the velocity channels from 
Vlsr = —7.5 km s~^ to 3.5 km s"^ on the blueshifted side, and from Vlsr =5.5 km s~^ to 
31.5 km s~^ on the redshifted side. Contour levels are 3a, 5a, and every 5a steps. The la 
levels are 1.4 Jy beam~^ km s~^ for the blueshifted component and 2.3 Jy beam~^ km s~^ for 
the redshifted component. The open circle is the field of view at 345 GHz, and the solid oval 
indicates the synthesized beam size. Note that the primary beam correction is not applied 
to the data. 
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Fig. 3. — A comparison of the CO 2-1 (a), CO 3-2 (b) and ^"^CO 2-1 (c) spectra observed 
with the SMA (black histograms) and the JCMT (grey histograms) at the pointing center. 
The SMA data were first primary beam corrected, and the spectra were smoothed to have 
the same angular resolution as the JCMT data, i.e., 14" for the CO 3-2 and 21" for the 
CO and ^^CO 2-1. The JCMT spectra are resampled to have the same velocity resolution 
as that of the SMA data and converted into brightness temperature by assuming the main 
beam efficienciies, //mb, to be 0.69 for the CO 2-1 and ^^CO 2-1, and 0.63 for the CO 3-2. 
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Fig. 4. — Velocity channel maps of the ^^CO 2-1 emission. The contours start at — 5cr and 
are drawn every 5a step until 30a, and every 10a in the range above 30cr. The la level is 
0.13 Jy beam~^. The crosses indicate the positions of the continuum sources A and B; the 
open circle is the SMA field of view at 230 GHz. Note that the primary beam correction is 
not applied to the data. 
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veloclty (km/s) 

Fig. 5. — Position-velocity (P-V) diagram of the CO 2-1 along a cut at P.A. = 90° through 
the position of source A. Contours are logarithmic at: 0.90, 1.43, 2.26, 3.58, 5.68, 9.00, 
14.3, 22.6, and 35.8 K. The horizontal thin dashed line is the position of source A; while the 
vertical thin dashed line labels the systemic velocity Vsys=4.0 km s~^. The green curve: wide- 
angle wind model curve produced by 2=30°, C=0.20 arcsec"^, and Vo=1.30 km s~^ arcsec"^. 
The black curve: wide-angle wind model curve produced by z=40°, C=0.10 arcsec~\ and 
fo=0.80 km s~^ arcsec"^. The filled square at the bottom left denotes the spacial and velocity 
resolution. 
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Fig. 6. — Outflow shell model curves overlaid onto the redshifted CO 2-1 emission. The 
green line and the black line delineate the curves with C = 0.20 arcsec"^ and C = 0.10 
arcsec"^, respectively. Both curves were produced by assuming the inclination angle of the 
axis from the plane of the sky to be 30°. 
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Fig. 7. — a: Curves of constant line ratios of CO 2-1/^^CO 2-1 (dashed lines) and constant 
CO 2-1 intensity (solid lines) as a function of the kinetic temperature and the H2 number 
density for dv/dr = 1000 km s"^ pc^^. The filled squares denote the corresponding observed 
values of bl and b2. b: Curves of constant line ratios of CO 3-2/CO 2-1 (dashed lines) 
and constant CO 2-1 intensity (solid lines). The locations of bl and b2 in this diagram 
correspond to the physical conditions derived from CO 2-1 and ^^CO 2-1 data. 
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Fig. 8.— Line profiles of tlie CO 2-1 (black line), ^^CO 2-1 (green line) and CO 3-2 (red 
line) at the compact components bl and b2. The position offsets of each component with 
respect to the phase tracking center are labeled in arcseconds at the top right. 



